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In Ref. equations of the modified geometrical optics for circularly polar- 



ized photon trajectories in a stationary spacetime are derived by using a (1+3)- 
decomposed form of Maxwell's equations. We derive the same results by using a 



four-dimensional covariant description. In our procedure, the null nature of the 
modified photon trajectory naturally appears and the energy flux is apparently null. 
We find that, in contrast to the standard geometrical optics, the inner product of 
the stationary Killing vector and the tangent null vector to the modified photon 
trajectory is no longer a conserved quantity along light paths. This quantity is fur- 
thermore different for left and right handed photon. A similar analysis is performed 
for gravitational waves and an additional factor of 2 appears in the modification due 
to the spin-2 nature of gravitational waves. 
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I. INTRODUCTION 

Light propagation in the gravitational field of a rotating body has been a topic of study 
in the past several years. One phenomenon of interest is rotation of the polarization vec- 
tor, known as gravitational Faraday rotation j^^- This effect does not manifest in the 
gravitational field of a non- rotating body, such as Schwarzschild spacetime joj , while it does 
occur for Kerr spacetime. This fact suggests that the existence of helicity-rotation coupling 
and the propagation of circularly polarized electro-magnetic waves depends on the helicity. 
Its occurrence has been confirmed by analyzing Maxwell's equations in curved spacetimes 



created by rotatin g bo dies |1CH12||. This effect is also discussed based on the gravitational 
Larmor's theorem 13|, uA . 

In Ref . [H , Frolov and Shoom reported that the spinoptics in a gravitational field created 
by a rotating massive compact object can be described by a modified geometrical optics 
approximation. They used a (H-3)-decomposed form of Maxwell's equations and also con- 
sidered a standard geometrical optics approach. In their setting of the base vector field for 
the circular polarization, a phase shift appears that depends on the helicity. They proposed 
a modification in which the ordering of the equations associated with the geometrical optics 
approximation is changed so that the phase shift is absorbed in the eikonal of the eikonal 
ansatz. This treatment also leads to a modification of the photon trajectory depending on 
the helicity. Using this procedure, scattering of circularly polarized light by a rotating black 



hole is discussed in Ref. 15 



In this paper, we derive the same equations of the modified geometrical optics as in 
Ref.Q by using another description in which four-dimensional covariance is maintained. In 
our procedure, we can easily see the four-dimensional picture of the photon propagation. In 
addition, this procedure can be easily applied to the case of gravitational waves, as will be 



explicitly shown (see also Ref. |14j|). 

This paper is organized as follows. In Sec. Ill At we review the standard geometrical optics 
approximation for Maxwell's equations in a stationary spacetime. We introduce a circular 
polarization base vector field in Sec. Ill Bl using an identical method to that in Ref. We 
then discuss the transport equation for the polarization vector using the circular polarization 
base vector field in Sec. Ill CI The photon trajectory, transport equation and energy fiux in 
the modified geometrical optics are given in Sec. IIIII In Sec. [IV]the procedure is apphed to 
the case of gravitational waves. Sec. |V]is devoted to a summary and discussion. 



II. STANDARD GEOMETRICAL OPTICS 



A. Trajectory, Transport Equation and Energy Flux 

In this paper we focus on a stationary spacetime manifold {A4,g), where is a four- 
dimensional manifold with a smooth Lorentzian metric g which has a smooth 1-parameter 
group G of the isometry generated by the Killing vector field ^. Following Ref. we write 
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the line element of the stationary spacetime as follows: 

ds^ = -h{dt - gidx'Y + h%dx'dx\ (1) 

where i, j run from 1 to 3 and h, cji and are functions of x*. The stationary Killing vector 
field is given by 

e = |, e^, = -h<Q. (2) 

Using the action of the isometry group G on A^, we can define the orbit space associated 
with the Killing vector field ^ as A/" := A^/G. We define the normalized Killing vector field 
u by 

:= ^7v^. (3) 
For later convenience, we define the projection tensor 7 by 

IpLv ■= Qi^v + Uf,u„. (4) 

Then, = h^ij and 7 gives the naturally induced metric on M . 

Note that in this paper, we consider the region in which h is positive definite. This con- 
dition may not be satisfied for regions within the ergosphere of a Kerr black hole. Therefore, 
as with the formalism in Refs. ll5| , our formalism cannot be straightforwardly applied to 



the ergoregion with the Killing vector field which is tangent to the world line of the static 
observer at the infinity. 

In Ref. Maxwell's equations are reduced to the master equations on the orbit space 
M with the metric 7. We do not follow the same procedure and instead use the four- 
dimensional covariant form of the equations. We consider the vector potential which 
satisfies the Lorenz gauge condition given by 

V^A'^ = (5) 

and the wave equation given by 

v^v.A^ - r;:a, = 0, (6) 



where R^y is the Ricci curvature tensor. Following the standard method (e.g. Refs. |l6L Il7]). 
we write the eikonal ansatz as follows: 

A, = {a, + eh^ + 0{e^)) e'"\ (7) 

where e is a book-keeping parameter that we take to be small during our manipulations; at 
the end of our calculations we reset it to e — )■ 1, so that S becomes the actual phase function. 

Substituting the ansatz ([7]) into Eq. (j5]), we obtain the following equation from the order 
of e-^: 

a^k^ = 0, (8) 

where kfj, is defined by 

k^ := V^S. (9) 
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From the order of e ^ in Eq. ([6]), we obtain 

k^kf, = 0. (10) 

Rewriting Eq. (|TOl) as 

n := ^g^-'V.SV.S = 0, (11) 

we can regard this equation as a Hamilton-Jacobi equation for S. Since the four-velocity of 
the corresponding dynamical system to Eq. ( ITTl) is given by k'^, we can regard the Hamilto- 
nian equation for the Hamiltonian ( ITTl) as the equation for the ray trajectory generated by 
k'^. The Hamiltonian equations are given by 

rv^A;'^ = 0. (12) 

Hence trajectories are given by null geodesies. This equation can be simply derived by 
differentiating Eq. (flO!) and using Eq. (|9]). 

The order of in Eq. ([6]) gives the following transport equation: 

k^V^a" + ^a^'V^k" = 0. (13) 

Following convention, we divide into the real amplitude a and the complex polarization 
vector as follows: 

= ai^ , i% = 1 , a G M, (14) 

where denotes the complex conjugate of i^. Then, contracting 1^ with the transport 
equation f|T3!) . from the real part, we obtain 

V^(a2p) = 0. (15) 

Substituting this equation into f|T3|) . we obtain 

rv^r = 0. (16) 

Eq. ( lT5l) describes the conservation of the photon number and Eq. ( IT6ll indicates that the 
polarization vector £^ is parallel-transported along the ray trajectory. 
The field strength F^j, of the vector potential ([7]) is given by 

F^, = Re {V^A, - VuA^} ~ 2aRe {2e'^fc[/,] } (17) 

at the leading order of the geometrical optics approximation, where square brackets denote 
anti-symmetrization. Then, the energy momentum tensor T^^ is given by 

T^"" = — [F^F"^ - -g^'^F^.F^'' 
47r V ^ 

= ^k^k" (1 - Re (e^^^^A^H) . (18) 

Averaging over several wavelengths, we obtain 

(T^^) = — fc^F. (19) 

This expression indicates that the energy flux is proportional to k^ and null at the leading 
order of the standard geometrical optics approximation. 
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B. Base Vector Fields 



Taking stationarity into account, we impose 

Cf:k^' = rV.A;^ - k'V.e = 0, (20) 

where is the Lie derivative with respect to ^. Using this equation and V^k^ = V^k^, we 
obtain 

v,(rM = o. (21) 

We define the frequency u as follows: 

to := -ek^,- (22) 
We introduce the spacelike unit vector along the ray direction n^, given by 

,= _ ^M. (23) 

u 

This satisfies 

n^'n^ = 1, nf'Uf, = 0. (24) 

To set an orthonormal base system, we define two additional unit spacelike vector fields 
and 63, given below. First, at a point, we set such that the following conditions are 
satisfied: 

gtive-A^^B = ^AB , u^e^j^ = n^e'^ = 0, (25) 

where A = 1, 2. Then, following Ref. we extend along the integral curve of by 
imposing the following condition: 

^ n'^D^e'X = -e\{n^Dxn^)n^' = n^'n^n^Dxe^, (26) 

where the action of on a vector field is defined by 

D,v, = l^lj'^pvx. (27) 

We can check that the condition Eq. ( l26l) is equivalent to Eq. (85) in Ref. which gives 
Fermi transport on (A/", 7). In addition, we extend the base vector fields along the integral 
curve of ^ by the Lie transport as follows: 

C^e^j^ = 0. (28) 

Then, Eq. fl25|) is satisfied at any point of the spacetime. 

Finally, we define the circular polarization base vector field as follows: 



= (e^ + iae^)/V2, (29) 

where cr = ±1 specifies circular polarization. Then, we have 

= V2{m^ + m^'), (30) 
= -iaV2{m'' -m''). (31) 

and also satisfy Eq. ( 126|) and are Lie transported along the ^ direction. 
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C. Parallel Transport of the Polarization Vector 

As shown in Eq. ( 1161) . the polarization vector is parallel-transported along the null 
geodesic generated by fc^. Using the circular polarization base vector field m'^, we can write 

t = m^'e''^, (32) 

where (p is a. real function of x*. Then, Eq. ( fT6l) can be rewritten as 

rV^(mV^) = ^ m^rV^(e^^) = -e^^FV^m^. (33) 

Contracting with m'^, we obtain 

ik^Vu^ = m^'k^Vufn^. (34) 

Using ( 123|1 . we find 

ijj CO 

OJ 00 

= — =171^12" D^rrin H m^^'^Vj^m^ 

Vh h 

"'m^m''V,^^, (35) 



h 

where we have used J-'„m^ = and C^m^ = 0. Since V^^jy is anti-symmetric, we obtain 

k'V,^ = a^eS^e^V^e^ 

= ^au.kxS^'-P^V^u^, (36) 

where e'^^'^^ is the completely anti-symmetric tensor with e'^^^s _ i / det g. Performing 
(1-1-3) decomposition, we can check that Eq. (136|1 is equivalent to Eq. (102) in Ref. [l[. This 
is the well-known gravitational analogue of the Faraday rotation |2|-|8|. 



III. MODIFIED GEOMETRICAL OPTICS 



A. Modification of the Eikonal 

The guiding principle of the modification is that the phase term e*''' should be included 
in the eikonal(see Eqs. ([7]), (JH]), and fl32|) ). that is, the modified eikonal S should be given 
by 

S S + if. (37) 
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Before modification, the Hamiltonian of the ray trajectory is given by 

n = ]^g^^k,K = \g^''V^SV,S. (38) 
The modification of the eikonal (|37|) and Eq. f p6|) suggest the following Hamiltonian: 

^ = ^^^'''(V^^ - av,)iy,S - av?^), (39) 

where we have defined 

:= ^e^,pxu''yPu\ (40) 

This expression was first derived by a group at Osaka City University 181] in a different way. 

Our aim is to modify the ordering of the field equations so that Eq. fl5^ is obtained. It 
will be seen in Eq. that our procedure eliminates the phase shift Eq. We do not 

change the form of the eikonal ansatz ([7]) but formally put the tilde " ~ " on all quantities, 
as follows: 

= + ehp + 0{^)) (41) 



To obtain the Hamiltonian (!39|) we change the orders of significance in the geometrical optics 
approximation by rewriting the gradient operator as follows: 

^- V/, - ie~Vv9^ + icrv9^. (42) 

Then, the Lorenz gauge equation and wave equation become 

\/pA^ = -> (V^ - %t-^a^^ + %o^^,)A^ = 0, (43) 
V^V^A^ = C(e°) ^ {V^-ie-^aip^ + iaip^){V'' -le-^aifi'' + iaip'')A^ = 0{e^). (44) 

This modification is trivial if we take e — )■ 1, but this enhances the effect of the circular 
polarization to the leading order. 

From the order of in Eq. fH5]) . we obtain 

a^q>. = 0, (45) 

where 

QtM = V^S^ - aLp^. (46) 

This means that the polarization vector a'^ is perpendicular to the ray direction given by q^. 
From the order of in Eq. (jH]), we have 

q% = 0. (47) 

This equation is identical to = and H is simply the Hamiltonian for the ray trajectory. 

In the same way as with Eq. ( l20l) . we extend the vector with the Lie transport along 
the integral curves of that is, 

C^q^" = 0. (48) 
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Then, we define the frequency uj as follows: 

^ := -eq,- (49) 

It should be noted that this frequency is not constant in general, in contrast to u = —^^k^. 
Since depends on the helicity cr, u also depends on a. We define the spacelike unit vector 
along the modified ray direction n'^ as follows: 

\/h 

W ■= X^q" - u^. (50) 

CO 

Following the same procedure as in Sec. Ill B[ we can define the modified circular polarization 
base vector m'^ associated with n^. 

From the order of in Eq. (HI]), we obtain 

g'V.S^ + ^a^V.q" + laq^'ip.W = 0. (51) 

We divide into the real scalar amplitude a and the circular polarization vector £^ := 
a'^/a = m^e*'^. Contracting with m , we obtain 

q'^V^a + iaq^Vy^p + am ^q'^V uTti^ + -^aVi^q^ + iaaq^ify = 0. (52) 
Similar to Eqs. ( 135|) and ( l36l) . the third term of this equation can be rewritten as 

rfi^q^V ufn^ = ^m^(n'' + M^) Vj^m^ 

yji 

-m,.u \ „m 



= -taqf^iff,. (53) 

Then, from the real and imaginary parts of Eq. ( 152|) . we obtain the following two equations: 

V^{a\n = 0, (54) 
g'^V^^ = 0. (55) 

Eq. fl5^ describes the photon number conservation and Eq. fl55l) indicates that the phase ^ 
is constant along the ray trajectory. Eq. (155!) is the desired result for the modification. 

From the Hamiltonian fl39l) . we obtain the following equation of motion for the ray tra- 
jectory: 

q'V.q'' = ar^q-", (56) 

where 

fi,u = ^fi^v - ^v^fi- (57) 
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Performing (H-3)-decomposition, we can derive Eq. (112) in Ref. [l|. 

We also perform the replacement P2|) in the expression f|T7j) . We obtain 

F^^ = Re { (V;, - ie'^oip^ + ioipf,) Ay - [V^ - ie'^oipu + A^} ~ 2aRe |^e*^g[^4] | 

(58) 

at the leading order of the modified geometrical optics approximation. Then, the energy 
momentum tensor T^^ is given by 

T^"" - ^q'^q'' (l - Re je^^'^^^^}) • (59) 
Averaging over several wavelengths, we obtain 

{Tn ^ fq'q". (60) 

This expression indicates that the energy flux is proportional to and null at the leading 
order of the modified geometrical optics approximation. 



IV. GRAVITATIONAL SPINOPTICS 

A. Standard Geometrical Optics 

The eikonal ansatz for the metric perturbation is given by follows: 

V = K- + eV + (^(e'))e*^/^ (61) 

Hereafter we work in the transverse-traceless gauge. The transverse gauge equation and the 
wave equation are given by 

V^/i^'^ = 0, (62) 
VpV^ V + '^Rp.xuh'^ = 0. (63) 

From the order of in Eq. (!62|) . we obtain 

a^uk" = 0, (64) 

where is given by ([9]). From the order of in Eq. ( l63l) . we obtain the same equation as 
Eq. (ITOl) and (|T2l) . From the order of in Eq. (1631) . we obtain 

k'Vpa^y + ia^.VpF = 0. (65) 

As in Eq. (1141) . we divide a^,y as follows: 

a^. = ai^, , = 1 , aeR. (66) 
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We can then obtain the following two equations: 

{a^k^) = 0, (67) 
k'V/^, = 0. (68) 

Eq. f l67|) describes the graviton number conservation and Eq. f l68|) indicates parallel transport 
of the polarization tensor along the null geodesic generated by k^. 

From Isaacson's formula 19|, |20|], the effective energy momentum tensor for gravitational 
waves can be written as 

^^^^""^^ = ^(Re{VAA}Re{V./.^^}>. (69) 
At the leading order of the geometrical optics approximation, we obtain 

This expression indicates that the energy flux of the gravitational waves is proportional to 
k^ and null at the leading order of the standard geometrical optics approximation. 



B. Base Setting and Parallel Transport of the Polarization Tensor 

Let us consider the base tensor fields for linear polarization tensors given by 

er = ^'5^''«, er = v^ei^e^, (71) 
where round brackets around indices denote symmetrization. These satisfy 

^Mi^^pAe^^'ex^ = 1 , ^^^^pAe7e+^ = 1 , ^^^^pAe^''e+^ = , gpxe'^^e^^ = e^l'e^ . (72) 

For the circular polarization specified by cr, we can define the polarization base tensor m^;^ 
by 

m'" = ^ (ef + ^aeT) • (73) 

This satisfies 

m^^m>"' = l, mf,^m>"' = l, fn^^m^"' = , m^'Plfip'' = iae^^e^ . (74) 
Using this circular polarization base tensor, we can write 

= m^,e'^. (75) 
Parallel transport of the polarization vector means that 

k^V p{m^,e''^') = ^ m^,fc^V,(e^'^) = -e'^^'k^V ^m^,. (76) 
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Contracting with m^'^, we obtain 

ik'^V^.tfj = m^^kPVpm^^. (77) 

Using Eq. f l2^ . we have 



00 ,,„ 00 



Vh ' h 

h' 

lau^kxeP^P^V^u^ = 2ik^V^ip, (78) 



where we have used e^rfV uTn^p = 0, C^mP" = and Eq. (!74|) . 



C. Modified Geometrical Optics 

Eq. f lTHj) suggests the Hamiltonian 

^gw = ^^7'''(V^5gw - 2(T^^)(V.5gw - 2a(^^) (79) 

for the modified geometrical optics. The extra factor 2 compared with Eq. ( I39ll is expected 
from the spin-2 nature of gravitational waves. We do not change the eikonal ansatz for the 
metric perturbation: 

hpu = {ap. + + 0(e2))e*^-/^ (80) 

As for Eqs. fH51) and flH|) . we rewrite the transverse gauge equation and the wave equation 
as follows: 

V^/i^" = ^ (V^ - 2ze-V^^ + 2ia^p)h^"' = 0, (81) 
VpVV = ^ (Vp - 2ie-V(/?p + 2icr¥?p)(V^ - 2te-^aipP + 2ia¥?'') ^ = C(e°)(82) 

From the order of in Eq. flHTj) . we obtain 

S/.^'P" = 0, (83) 

where 

Pt, = V^S'gw - 2aLp^. (84) 
From the order of in Eq. fl82p . we obtain 

P% = 0. (85) 
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This equation is identical to "Hg^ = 0. In the same way as for IIII At we impose 

C^p'' = (86) 

and define the frequency Wg-w as follows: 



gw •- 



(87) 



Then, we define the spacelike unit vector along the modified ray direction n^^ by 



gw ~ 



-pf" - u^. (88) 



We can then obtain the modified circular polarization base tensor m^y associated with n^. 
From the order of in Eq. ( 182|) . we obtain 

P''^ lAp^v + \^i,v'^ pV^ + 2iapPippd^y = 0. (89) 

We replace a^^ by am^ye^'^ . Contracting Eq. ( I89l) with fh^ , we obtain 

p^'Wya + iap'^Vutjj + am^^pPV pm>"' + ]^aV ^p"" + 2iaap''ify = 0. (90) 



Similar to Eq. (!78i) . the third term of this equation can be rewritten as 



m p'^Wpm^y = ---j=m^y{nl^ + u'')Vpm 



'm^yU''Vpm 



= -2tap''ip^. (91) 

Then, from the real and imaginary parts of Eq. f l90p . we obtain the following two equations: 

V^(«y) = 0, (92) 
P^V^tP = 0. (93) 

Eq. (1921) describes the graviton number conservation and Eq. (1931) indicates that the phase 
is constant along the ray trajectory. The equation of motion for the ray trajectory can be 
derived from the Hamiltonian (1791) as follows: 

p-'V.pP = 2ary, (94) 

where is defined in Eq. (|57|) . 

We also perform the replacement (142|) in the expression (l69l) . and find that 

{T^') - ^"VV. (95) 

This expression indicates that the energy fiux is proportional to p'^ and null at the leading 
order of the modified geometrical optics approximation. 
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V. SUMMARY AND DISCUSSION 

Using a four- dimensional covariant description, we have derived the equations of the 
modified geometrical optics in a stationary spacetime, previously derived in Ref . jl| . In the 
modified geometrical optics, the three-dimensional photon trajectory is modified depending 
on the photon helicity. In Ref. [l|, the authors used a reduced form of Maxwell's equations on 
the three-dimensional orbit space associated with the stationary Killing vector field. In this 
description, the four-dimensional picture is not clear. In contrast, in our procedure, the null 
nature of the photon trajectory in the modified geometrical optics naturally appears and the 
energy flux is apparently null. We can also see that, in contrast to the standard geometrical 
optics, the inner product of the stationary Killing vector and the tangent null vector to the 
modified photon trajectory is no longer a conserved quantity along light paths. This quantity 
is furthermore different for left and right handed photon. The same procedure can be easily 
applied to the case of gravitational waves and we found that an additional factor of 2 appears 
in the modification between the circularly polarized photon and the graviton because of the 
spin 1 and 2 nature of electro-magnetic waves and gravitational waves, respectively. 

It is clear that the origin of the modification is in the choice of the circular polarization 
base vector field. Following Ref. we have taken a circular polarization base vector field 
based on the Fermi transport along the photon trajectory projected on the three-dimensional 
orbit space. However, it is still not clear whether this choice of base vector field is valid 
for describing the circularly polarized photon trajectory. Other choices of the base vector 
field may give rise to different modifications. Clearly we need to justify the choice of base 
vector field based on observations. This could be done by comparing the electro-magnetic 
fields given by the modified geometrical optics with those given by directly solving the wave 
equations. We leave this issue to a future work. 
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